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^Studies  were  directed  at  direct  measurement  of  the  acoustic  admittance  of 
burning  solid  propellant  with  concurent  pressure  and  optical  anemometry 
measurements,  internal  ballistics  of  spin  stabilized  nozzleless  rocket  motors, 
and  the  application  of  Hilbert  transforms  to  the  computation  of  the  imaginary 
part  of  the  pressure  coupled  response  from  data  for  the  real  part.  An  optical 
anemometer  based  on  an  Argon  ion  laser  and  tracker  electronics  and  a  servo- 
positioned  window  strand  burner  based  on  a  Princeton  window  burner  were 
constructed.  The  servo-positioned  strand  burner  operated  successfully  and  \ 
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Doppler  burst  phenomena  Mere  observed.  However,  tracker  lock-on  was  not 
achieved.  The  spin  stabilized  nozzleleas  analysis  showed  that  gas  dynamic  spin 
effects  (as  contrasted  with  acceleration  induced  burning  rate  changes)  were 
small  for  tangential  Mach  numbers  below  0.1.  The  Hilbert  transform  study  showed 
that  the  imaginary  part  of  the  response  function  could  be  computed  from  the  real 
part  with  acceptable  accuracy  for  frequencies  below  one-third  of  the  maximum 
frequency  of  the  data.  . 
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SECTION  1 
INTRODUCTION 


The  design  process  employed  with  solid  propellant  rockets  involves  the 
analysis  of  quasi-steady  performance,  structural  integrity,  combustion 
stability,  and  cost  for  candidate  designs.  Performance  and  stability  are  ana¬ 
lyzed  with  codes  such  as  SPP^-  and  SSP^ .  These  codes  require  data  on  pro¬ 
pellant  properties  and  the  motor  geometry.  Typically,  the  performance  code 
(SPP)  supplies  the  stability  code  (SSP)  with  the  required  grain  geometry,  the 
local  burning  rate,  and  the  mean  port  flow.  Therefore,  if  appropriate  pro¬ 
pellant  parameters  are  available,  these  codes  can  evaluate  both  the  perform¬ 
ance  (pressure-  and  thrust-time  histories)  and  the  propensity  for  instability 
(modal  growth  constants)  of  candidate  designs. 

The  experimental  determination  of  propellant  parameters  for  the  stability 
analysis  is  difficult,  uncertain  relative  to  velocity-coupled  and  distributed 
combustion  driving,  and  expensive.^  This  situation  and  the  need  for  these 
design  data  have  led  to  a  search  for  and  development  of  improved  diagnostics. 
Two  direct  measurement  techniques  have  been  examined:  microwave  inter¬ 
ferometry^  and  optical  anemometry.^  The  former  seeks  to  measure  the 
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thickness  of  the  solid  propellant  instantaneously  so  that  instantaneous 
burning  rate  and  subsequently  the  response  function  can  be  computed  from  its 
temporal  derivative.  The  latter  seeks  to  measure  the  instantaneous  gas  velo¬ 
city  at  the  outer  edge  of  the  gas  phase  reaction  zone,  so  that  the  acoustic 
admittance  can  be  determined  directly.  From  continuity  the  instantaneous  gas 
velocity  is  related  to  the  instantaneous  burning  rate  by* 


ug  ~  Pcrb/pg  • 


(1) 


Since  the  order  of  magnitude  of  pc/pg  is  roughly  100,  gas  velocity  fluc¬ 
tuations  will  be  roughly  two  orders  of  magnitude  larger  than  burning  rate 
fluctuations.  Moreover,  thickness  changes  aL  that  must  be  measured  by  the 
interferometer  are  of  magnitude  rbAt,  where  rb  is  amplitude  of  the  burning 
rate  fluctuation,  and  At  is  the  period  of  the  oscillation.  Now 
9(3  ~  Rpfjj  (p/p).  Thus,  the  thickness  change  due  to  nonsteady  burning  in  one 
cycle  is 

AL  ~  Rprb(p/p)/f.  (2) 


From  (1)  and  (2) 


A  A 

Ug  ~  RpFb(p/p)  PC/pg. 


(3) 


*  Continuity  gives  3p/3t  +  v*pu  =  0.  Making  gas  velocity  measurements  near 
planes  of  symmetry  reduces  7*pCf  to  dpu/3x;  v,  w,  3pv/3y  and  3pw/3z  are  small 
becauae  of  symmetry.  Consequently,  continuity  becomes  3p/3t  +  3pu/Sx  ~  0. 
Integrating  and  applying  Liebnitz  rule  and  the  mean  value  theorem  for 
integrals  gives  Aydp/dt  +  uQpa  -  rbpe  s  0.  For  small  oscillations  and  Ax  ~  0* 
(1)  results. 
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Comparing  (2)  and  (3)  it  i8  clear  that  the  magnitude  of  the  length  pertur¬ 
bation  that  must  be  measured  (aL)  is  at  least  two  orders  smaller  than  the 
magnitude  of  the  velocity  perturbation  (ug);  at  1000  Hz  it  is  five  orders  of 
magnitude  smaller.  Moreover,  the  magnitude  of  AL  decreases  as  frequency 
increases.  For  typical  values  (Rp  ~  4,  r^  ~  1  cm/s,  p/p  ~  0.1,  f  ~  1  kHz, 
pc/pg  ~  102)  aL  ~  4  x  10~4  cm  =  4ym  and  Ug  ~  40  cm/s  (the  mean  gas  velocity  is 
100  cm/s).  For  10  gH  microwave  signals  the  wavelength  is  roughly  3  cm.  Thus, 
measurement  of  response  functions  at  1  kHz  requires  a  phase  measurement 
precision  of  the  order  of  one  part  in  10^  to  reach  one  percent  accuracy.  This 
is  at  the  limit  of  current  technology.  Consequently,  microwave  measurements 
in  the  10  kHz  frequency  regime  are  doubtful.  By  contrast,  velocity 
fluctuations  on  the  order  of  40  cm/s  about  a  mean  velocity  of  100  cm/s  to 
frequencies  on  the  order  of  10  kHz  are  well  within  the  state  of  the  art  of 
optical  anemometry.  Accordingly,  anemometry  represents  perhaps  the  best 
approach  to  high-frequency  response  function  measurement. 

The  above  discussion  pertains  to  measurements  that  characterize  the 
coupling  of  energy  release  to  pressure  oscillations  in  the  first  one  milli¬ 
meter  or  so  from  the  burning  surface.  However,  with  some  propellants  (e.g., 
nitr amines)  the  gas  phase  reactions  are  slower,  so  that  the  energy  release  is 
spread  out  over  an  appreciably  longer  distance  (10*  -  10*  times  larger).  In 
these  cases,  distributed  combustion  driving  must  be  considered.  These 
situations  are  of  current  and  future  interest  because  nitramines  will  be  the 
basis  for  min-smoke  composite  propellants.  In  these  situations  one  also  wants 
to  know  the  way  the  quasi-steady  energy  release  is  distributed,  so  that  it  can 
be  correlated  with  that  fraction  of  the  chemical  energy  release  that  ia 
coupled  to  pressure  oscillations.  This  correlation  would  permit  quasi-steady 
calculations  (available  technology  for  CP  grains**)  to  form  the  basis  for  com¬ 
puting  distributed  combustion  driving. 


6.  Beddini,  R.  A.,  "Reacting  Turbulent  Boundary  Layer  Approach  to  Solid 
Propellant  Erosive  Burning,"  AIAA  Journal,  16,  No.  9,  Sept.  1978,  pp.  898-905. 
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For  quasi-steady  flow  along  planes  of  symmetry,  transverse  variations  and 
velocity  components  are  small.  Thus  continuity  leads  to  (1),  which  can  be 
rewritten  using  the  perfect  gas  law  as 


RgTg  -  P  Ug/r^pQ  •  (4) 

Clearly,  with  concurrent  pressure,  gas  velocity,  and  burning  rate  measurements 
RgTg  can  be  determined.  In  the  absence  of  transverse  gradients,  the  quasi¬ 
steady  energy  equation  becomes 


dH  =  c  dT  +  dQ„  +  du„2/2  a  0  . 
p  g  R  g 


Thus,  the  local  energy  release  is 


^R 


-c  dT  -  duJ/2 

p  g  g 


3_ 

3x 


(5) 


(6) 


For  ideal  gases  cp  =  Rg  k/(k-l).  Moreover,  in  the  distributed  energy  release 
region  k  does  not  vary  atrongly.  Therefore,  with  (4) 


dQR  “ 


klT  d(Pug)/rbpc  •  dug2/Z  ‘ 


(7) 


Once  again,  it  i8  clear  that  local  qua8i-steady  heat  release  can  be  measured 
with  concurrent  burning  rate,  pressure,  and  gas  velocity  measurements  along 
the  symmetrical  stream  tube.  The  variation  of  specific  impulse  and 
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characteristic  velocity  with  characteristic  length  can  be  estimated^  with 
RgTg  known  as  a  function  of  distance  (and  time*)  from  the  burning  surface. 

The  mechanical  energy  equation^  shows  that  the  conversion  of  internal 
energy  into  mechanical  energy  is 


dME/dt  s  p  $*u  . 


(8) 


Therefore,  for  a  stream  tube  along  planes  with  transverse  symmetry  the  net 
mechanical  energy  added  to  the  nonsteady  field  is 


d(ME)'/dt  =  p’3u'/3x  •  (9) 

For  a  bulk  mode  oscillation  as  might  occur  in  a  strand  burner  this  is 

d(ME) '/dt  =  3(p*u' )/3x  *  £p'u'(x+Ax)  -  p'u'Cxy  J/ax  .(10) 

Thus  distributed  combustion  driving  can  be  measured  by  mapping  out  p'u'  as  a 
function  of  distance  from  the  burning  surface.  This  can  be  accomplished  with 
simultaneous  pressure  and  velocity  measurements  at  either  one  location  that  is 


7.  Williams,  F.  A.,  Barrere,  M. ,  and  Huang,  N.  C.,  Fundamental  A8oect8  of 
Solid  Propellant  Rocketa.  AGAROograph  116,  October  1969,  pp.  173-177. 

8.  Bird,  R.  B.,  Stewart,  W.  E.,  and  Lightfoal,  E.  N.,  Transport  Phenomena. 
(John  Wiley  and  Sons,  New  York,  1960)  pp.  81-82. 


*  time  =  Jx  u"1 
o  9 
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slowly  traversed  relative  to  the  burning  surface*  or  at  two  closely  spaced 
locations  that  are  slowly  traversed  relative  to  the  burning  8urface.* 
Concurrent  pressure  and  velocity  measurements  may  also  be  employed  to  deter¬ 
mine  the  distributed  combustion  driving. 

In  summary,  it  is  seen  that  concurrent  pressure,  burning  rate,  and  gas 
velocity  measurements  in  the  core  region  of  a  strand  plume  can  be  employed  to 
determine  the  following: 

•  acoustic  admittance  of  the  burning  surface, 

•  local  temperature, 

9  local  energy  release, 

•  local  distributed  combustion  driving,  and 

•  specific  impulse  and  characteristic  velocity  and  their 
variation  with  residence  time. 

The  fact  that  these  measurements  can  be  made  concurrently,  with  reasonable 
precision,  and  to  frequencies  in  the  10  kHz  domain  indicates  that  the  optical 
anemometry  approach  offers  appreciable  diagnostic  advantages  relative  to  other 
methods . 

Caveny,  Collins,  and  Cheng^  showed  that  direct  measurements  of  the 
acoustic  admittance  of  a  burning  propellant  surface  are  feasible  with  optical 
anemometry.  However,  these  measurements  were  performed  with  reasonably  short 
strands  (~  5  cm).  To  achieve  increased  accuracy  and  the  potential  for 
acquiring  other  data,  a  longer  burn  time  is  required.  This  program  was 
primarily  aimed  at  repeating  the  Caveny,  et.  al.  experiments  with  a  longer 
burn  time. 


*  for  8inusoidal  oscillations  this  might  be  accomplished  with  a  lock-on 
amplifier  (p *  employed  as  reference)  having  an  appropriate  time  constant. 


SECTION  2 
ACCOMPLISHMENTS 


Progress  has  been  made  in  the  following  areas: 

•  construction  and  demonstration  of  a  servopositioned  strand  burner 
capable  of  burning  1/4  x  1/4  x  8  in.  strands, 

•  construction  of  an  optical  anemometer  using  an  argon  ion  laser  and 
tracker  electronics, 

o  analysis  of  the  effects  of  spin  on  the  idealized  performance  of 
nozzleless  rocket  motors, 

•  interpretation  of  spontaneous  transitions  to  nonlinear  instability  in 
forced  longitudinal  wave  motor  tests,  and 

«  application  of  Hilbert  transforms  to  the  computation  of  the  imaginary 
part  of  the  pressure-coupled  response  function  from  the  real  part  in 
general  situations. 

The  first  two  items  are  the  major  component  parts  of  apparatus  for  the  direct 
measurement  of  acoustic  admittance  of  burning  solid  propellant. 

Servopositioned  Strand  Burner 

A  standard  optical  strand  burner,  digital  stepping  motor,  and  electronic 
control  system  were  obtained  from  Princeton  University.  These  items  were 
incorporated  into  a  servopositioned  strand  burner  capable  of  burning  1/4  x  1/4 
in.  strands  up  to  eight  inches  long.  The  design  of  this  apparatus  is 
described  in  a  technical  note  submitted  for  publication  to  The  Journal  of 
Spacecraft  and  Rockets. 

Optical  Anemometer 

A  single-channel  Laser  Doppler  Velocimeter  (LDV)  system  was  assembled. 
Most  of  this  equipment  was  furnished  by  the  University  of  Dayton  Research 
Institute  (UDRI)  for  this  program.  The  LDV  system  includes  a  CW  argon  ion 
laser  and  a  TSI  1090-1A  tracker.  Successful  trials  of  this  system  have  been 
run  in  the  burner  geometry  (including  windows)  using  both  a  calibration  wheel 
and  smoke  88  a  target.  Initial  checks  on  propellant  flames  suffered  from 
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insufficient  signal-to-noise  ratio  for  useful  results.  It  is  believed  that 
the  system  has  now  been  improved  to  the  point  where  success  with  propellant 
strands  may  be  anticipated. 

Schlieren  and  Shadow  Photography  of  Burning  Rocket  Propellant 

Several  attempts  were  made  to  record  schlieren  photographs  and  sha¬ 
dowgraphs  of  the  gases  evolved  by  a  6-mm  square  by  50-mm  long  section  of 
rocket  propellant  ignited  at  one  end.  A  pulsed  metal  vapor  laser  producing 
20-ns  wide  pulses  of  light  was  used  as  an  illuminating  source.  The  laser  was 
fired  by  a  triggering  system  arranged  so  that  it  fired  synchronously  with  the 
motion  of  the  frames  on  a  16-mm  Fastax  camera.  Kodak  RAR  2496  film  in 
125-foot  rolls  was  used  as  the  recording  medium.  The  Fastax  camera  was 
capable  of  reaching  framing  rates  of  8  kHz  near  the  end  of  the  roll. 

The  light  from  the  laser  was  focused  onto  a  pinhole,  collimated  by  a 
lens,  and  passed  around  the  sample.  The  light  was  condensed  by  a  second  lens 
onto  the  film.  The  object  to  image  size  ratio  was  approximately  3:1.  The 
maximum  film  velocity  in  the  camera  was  60  m/s.  Since  the  exposure  time  was 
only  20  ns,  no  motion-compensation  prism  was  used  in  the  camera  because,  with 
such  a  short  exposure  time,  the  maximum  blur  of  the  image  was  only  1.2  urn,  far 
below  the  resolution  limit  of  any  16-mm  film.  The  sample  was  ignited  and  the 
camera  started  a  short  time  later,  because  the  rate  at  which  the  fuel  was  con¬ 
sumed  was  very  slow  compared  to  the  framing  rate  of  the  camera. 

In  general,  shadowgraphs  are  capable  of  defining  the  boundary  between  two 
volumes  of  gas  having  different  indices  of  refraction,  but  they  fail  to  yield 
quantitative  information  because  the  fringes  observed  depend  upon  the  second 
derivative  of  this  index.  They  are  useful,  however,  because  of  their  sharp 
definition  of  boundaries.  Schlieren  photographs  are  more  quantitative  because 
their  fringes  depend  upon  the  integrated  optical  path  length  through  the 
sample.  Regardless  of  the  fact  that  the  beam  stop  in  a  schlieren  system  elim¬ 
inates  a  considerable  amount  of  the  light  passing  around  and  through  the 
sample,  both  types  of  photographs  can  be  recorded  with  the  same  amount  of 
laser  light. 

Successful  schlieren  and  shadow  movies  were  recorded,  despite  the 
brilliance  of  the  flame  produced  by  the  combusting  fuel.  It  was  found  that  it 
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was  not  necessary  to  eliminate  the  yellow  light  produced  by  the  fuel  (although 
this  could  have  been  accomplished  easily  by  the  use  of  an  optical  filter) 
because  the  brightness  of  the  laser  was  so  much  greater  than  that  of  the 
flame.  In  general,  the  quality  of  the  films  produced  was  poor.  The  volume  of 
gas  produced  by  the  fuel  was  very  large  and  it  spread  far  beyond  the  three- 
inch  diameter  viewing  area  that  was  available.  The  fuel  burned  slowly  on  the 
time  scale  of  the  photographs.  Therefore,  future  attempts  to  record  the 
evolved  gases  on  a  global  scale  should  be  performed  at  a  lower  framing  rate, 
perhaps  as  slow  as  1  kHz.  Two  experimental  difficulties  led  to  the  poor 
quality  of  the  photographs  recorded:  (1)  in  the  case  of  the  shadowgraphs,  the 
system  sensitivity  was  simply  too  high  and  (2)  for  the  process  of  recording 
the  schlieren  photographs,  a  satisfactory  beam  stop  was  not  available. 

In  subsequent  work  on  other  combusting  systems,  these  deficiencies  were 
rectified  but  could  not  be  applied  to  the  photography  of  the  solid  rocket 
propellant.  The  sensitivity  of  a  shadowgraph  depends  on  the  percentage  of  the 
light  deflected  by  the  subject  which  is  allowed  to  reach  the  film  and  upon  the 
degree  to  which  the  subject  is  in  focus.  Because  of  the  large  size  of  the  gas 
cloud  from  the  fuel,  it  is  impossible  to  have  the  entire  cloud  in  focus  at  one 
time,  so  the  sensitivity  of  the  shadowgraphs  must  be  reduced  by  the  insertion 
of  a  stop  in  the  focal  plane  of  the  camera  lens  to  limit  the  amount  of  off- 
axis  light  reaching  the  film.  Such  an  aperture  was  not  available  when  the 
fuel  photographs  were  taken.  The  problem  associated  with  the  schlieren  pho¬ 
tography  was  created  by  the  use  of  coherent  laser  light  whose  diffraction  and 
interference  properties  at  a  sharp  beam  stop  edge  are  different  from  those  of 
incoherent  light.  The  combusting  rocket  fuel  was  the  first  example  of  the 
severe  gradients  present  in  a  harsh  flame. 

The  experience  gained  from  the  attempts  to  photograph  this  system  was 
useful  in  subsequent  schlieren  and  shadow  photography  performed  on  two  com¬ 
bustors  in  the  Aero-Propulsion  Laboratory  at  Wright-Patterson  Air  Force  Base. 

Spin  Stabilized  Nozzleless  Rocket  Motors 

An  integral  method  analysis  was  completed  for  the  effects  of  spin  on  the 
idealized  performance  of  nozzleless  rocket  motors.  Results  show  negligible 
effects  for  tangential  Mach  numbers  (tangential  velocity  of  burning 
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surf ace/ acoustic  speed)  below  0.1.  An  article  describing  this  work  has  been 
submitted  for  publication  to  the  forthcoming  AIAA  Joint  Propulsion  Conference. 


Spontaneous  Transition  to  Nonlinear  Instability  in  Forced  Longitudinal  Wave 
Motor  Experiments 

Literature  concerning  transition  to  nonlinear  instability,  acoustic 
erosivity,  erosive  burning,  and  cold  flow  studies  was  combined  with  FLW  motor 
results  to  explain  the  cross  sectional  geometry  dependence  of  the  transitional 
phenomena.  The  results  of  this  "inter-group  study"  were  reported  to  the 
AIAA/SAE/ASME  Joint  Propulsion  Conference  in  July  1981. 

Application  of  Hilbert  Transforms  to  Response  Function  Data 

Pressure-coupled  response  function  data  obtained  from  T-  and  rotating 
valve  burners  are  constrained  to  the  real  part.  However,  knowledge  of  the 
sign  of  the  imaginary  part  is  very  useful  for  design  purposes;  it  can  be 
employed  to  determine  the  sign  of  the  velocity-coupled  driving  contribution. 
However,  no  oeneral  methodology  was  known  to  predict  the  imaginary  part  from 
the  real  part.  The  applicability  of  Hilbert  transform  relations  to  this  task 
was  explored.  The  results  of  this  study  were  submitted  for  publicstion  to  the 
forthcoming  AIAA  Joint  Propulsion  Conference. 


10 


4 


PUBLICATIONS 


Glick,  R.  L.,  Micci,  M.  M.  and  Caveny,  L.  H.,  Transition  to  Nonlinear 
Instability  in  Solid  Propellant  Rocket  Motors.  AIAA-81-1520. 

Glick,  R.  L.,  Lincks,  D.  R.,  and  Yeruahalmi,  S. ,  "Servopositioned  Strand 
Burner,"  submitted  to  The  Journal  of  Spacecraft  and  Rockets  for 
publication,  January  1982. 

Yerushalmi,  S.  and  Glick,  R.  L.,  "On  the  Internal  Ballistics  of  Spin 
Stabilized  Nozzleless  Rocket  Motors,"  submitted  for  publication  to  the 
1982  AIAA  Joint  Propulsion  Conference,  December  1981. 

Becker,  R.  J.  and  Glick,  R.  1.,  "Application  of  Hilbert  Transforms  to 
Response  Function  Data,"  submitted  for  publication  to  the  1982  AIAA  Joint 
Propulsion  Conference,  December  1981. 


11 


NOMENCLATURE 


Cp  specific  heat  at  constant  pressure 

f  frequency 

dH  total  enthalpy  change 

k  specific  heat  ratio 

AL  strand  length  change  during  one  oscillation  due  to 

nonsteady  burning 

ME  mechanical  energy 

p  pressure 

dQp  heat  release  from  chemical  reaction 

rj,  burning  rate 

Rg  gas  constant 

Rp  pressure-coupled  response  function 

t  time 

Tg  gas  temperature 

Ug  gas  velocity  normal  to  burning  surface 

u  gas  velocity 

X  distance  from  burning  surface 

pc  propellant  density 

pg  gas  density 

(  )'  denotes  a  fluctuating  term 

(  )  denotes  a  mean 

A 

(  )  denotes  an  amplitude 
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